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Detecting dark energy dynamics is the main quest of current dark energy research. Addressing 
the issue demands a fully consistent analysis of CMB, large scale structure and SN-Ia data with 
multi-parameter freedom valid for all redshifts. Here we undertake a ten parameter analysis of 
general dark energy confronted with the first year WMAP, 2dF galaxy survey and latest SN-Ia 
data. Despite the huge freedom in dark energy dynamics there are no new degeneracies with 
standard cosmic parameters apart from a mild degeneracy between reionisation and the redshift of 
acceleration, both of which effectively suppress small scale power. Breaking this degeneracy will 
help significantly in detecting dynamics, if it exists. Our best-fit model to the data has significant 
late-time evolution at z < 1.5. Phantom models are also considered and we find that the best-fit 
crosses w = — 1 which, if confirmed, would be a clear signal for radically new physics. Treatment 
of such rapidly varying models requires careful integration of the dark energy density usually not 
implemented in standard codes, leading to crucial errors of up to 5%. Nevertheless cosmic variance 
means that standard ACDM models are still a very good fit to the data and evidence for dynamics 
is currently very weak. Independent tests of reionisation or the epoch of acceleration (e.g. ISW- 
LSS correlations) or reduction of cosmic variance at large scales (e.g. cluster polarisation at high 
redshift) may prove key in the hunt for dynamics. 

PACS numbers: 



I. INTRODUCTION 

Cosmological observations suggest that the Universe 
is dominated by an exotic form of matter which is re- 
sponsible for the present phase of accelerated expansion 
[E H> IE 111 ' Several scenarios have been proposed to ac- 
count for the observations, but the nature of this dark en- 
ergy still remains unknown. The simplest minimal model 
to fit the experimental data assumes the presence of a 
cosmological constant A, representing the vacuum energy 
contribution to the spatial curvature of the space-time. 
In spite of the success of this concordance ACDM model, 
there is no reasonable explanation why the observed value 
of A is extremely small compared to particle physics ex- 
pectations p. 

Alternatively, a light scalar field, called quintessence, 
rolling down a fiat effective potential has been proposed 
to account for the missing energy in the Universe |a, ■ 
In particular quintessence models manifesting 'tracker' 
properties allow the scalar field to dominate the present 
Universe independently of the initial conditions [a, |3- 
Different realizations of the original quintessence idea 
have been studied in the literature including the possibil- 
ity of a scalar field evolution driven by a non-canonical 
kinetic term 10] and a non-minimal coupling between 
quintessence and dark matter [ll|, |l3, EE, 13 . On the 
other hand unified models of dark matter and dark en- 
ergy have been considered [13, 0, C3 • Despite the pro- 
liferation of scalar field models of dark energy, we still 
lack of a fully consistent particle physics formulation. 



Nonetheless distinguishing between a dynamical form of 
dark energy and a cosmological constant is of immense 
importance as it would give us a hint on the nature of 
this component. The recent WMAP satellite measure- 
ments of the Cosmic Microwave Background, by provid- 
ing an accurate determination of the anisotropy power 
spectrum, offer the opportunity to have a better insight 
into the physics of the dark energy. The quintessence 
hypothesis has been tested with different methods using 
pre- WMAP CMB data and SN-Ia data or the 2dF power 
spectrum |ll HUll El IH E^. These analysis have 
constrained the dark energy equation of state w with- 
out ruling out the possibility of a time dependence. In 
this article we carry out an analysis of the time evolu- 
tion of the dark energy equation of state. Our aim is 
to constrain a set of parameters characterizing the dark 
energy properties and the standard cosmological param- 
eters by perfor ming a likelihood analysis of the WMAP 
first year data ^| and the SN-Ia luminosity distance 
measurements llU25J. The paper is organized as follows: 
in Section Inland the appendices we describe the method 
and the data, in Section IlIII we explain the evaluation of 
the likelihood, in Section}^ we discuss our results and 
finally in Section we present our conclusions. 



II. METHOD AND DATA 

The most general way of constraining the time evolu- 
tion of the dark energy equation of state would require 



us to consider a completely free function w{z). As this 
corresponds to an infinite number of new degrees of free- 
dom, we have to simplify the problem. We use a physi- 
cally motivated parametrisation where w — p/ p is defined 
by its present value, wq, its value at high redshift, Wm, 
the value of the scale factor where w changes between 
these two values, at and the width of the transition, A. 
Namely: 



w{a) ^ Wo + (wm - wo)r(a, at, A) 



(1) 



where F, the transition function, has the limits T(a — 
0) = 1 and r(a = 1) = and varies smoothly be- 
tween these two limits in a way that depends on the 
two parameters at and A (see figure^. Such a choice 
has been shown to allow adequate treatment of generic 
quintessence and to avoid the biasing problems inherent 
in assuming that w is constant. Two choices for F have 
been given in the literature [23, llg as discussed in Ap- 
pendix^ Here we use the form advocated in J2y|. 




FIG. 1: Schematic plot of the equation of state paratrisation 
Eq. dl. 

Using this general prescription has a profound advan- 
tage in attempts to detect dark energy dynamics since, 
unlike simpler parametrisations based on only one or 
two variables, it can accurately describe both slowly and 
rapidly varying equation of states [23. Detecting dark 
energy dynamics and distinguishing it from a cosmolog- 
ical constant is difficult and is clear only when there are 
rapid, late-time, changes in w 28j, which then needs a 
formalism capable of describing such rapid transitions. 

In order to compute the CMB power spectra, we use a 
modified version of the CMBfast Boltzmann solver [23| . 

This is a non-trivial step in the case where F changes 
rapidly (such as in our best- fit model!). In fact using a 
numerical method that is not able to track rapid transi- 
tions can lead to errors significantly larger than the error 
bars on the data, of order 5%, and consequently lead 
to completely wrong results. Our tests are described in 
detail in Appendix IbI and lO 

Several degeneracies amongst the cosmological param- 
eters prevent us from accurately constraining cosmologi- 
cal models using CMB data only. Specific features of the 



anisotropy power spectrum can provide information on 
particular combinations of the cosmological parameters. 
For instance the relative height of the Doppler peaks de- 
pends on the baryon density and the scalar spectral in- 
dex. In order to break such degeneracies it is necessary to 
add external information. Since our goal is to constrain 
the properties of the dark energy, in a flat geometry the 
main limitation comes from the geometric degeneracy be- 
tween wq, the dark energy density ftoE and the Hubble 
parameter h. This degeneracy can be broken by assum- 
ing an HST prior on the value of h [30j or/and combining 
the CMB with other data sets such as the matter power 
spectrum measurements from the 2dF galaxy survey |3l| 
or the SN-Ia data. In our analysis we use the "gold" sub- 
set of the recent compilation of supernova data of [23 
in addition to the WMAP TT and TE spectra. An im- 
portant point which we want to stress here is that CMB 
and SN data can be treated at a fundamental level with- 
out any prior assumption on the underlying cosmological 
model. For instance, this is the case for the matter power 
spectrum data from galaxy surveys which implicitely as- 
sumes a ACDM model when passing from redshift space 
to real space. For this reason we add the 2dFGRS large 
scale structure data only in order to check the stability 
of our results. We also remark that the use of secondary 
observables such as the age of the Universe, the size of 
the sound horizon at the decoupling, the clustering am- 
plitude cTg (as quoted by the WMAP-team) or the growth 
factor of matter density perturbations should not be used 
without thought to infer constraints on the dark energy 
since their quoted value is usually derived by implicitly 
assuming a ACDM cosmology. This leads to biased re- 
sults since these observables depend on the nature of the 
dark energy [23, 113, 113 ■ In principle CMB constraints 
can be easily added by using the position of the Doppler 
peaks, which provide an estimate of the angular diameter 
distance to the last scattering surface. But it is a well 
known fact that pre-recombination effects can shift the 
peaks from their true geometrical position |32l | . 

The Integrated Sachs- Wolfe (ISW) effect also induces 
an additional shift in the position of the first peak, in a 
way that is strongly dependent on the evolution of the 
dark energy equation of state and is generally larger than 
in ACDM models |23|. Therefore a consistent dark en- 
ergy data analysis of the CMB indeed requires the com- 
putation of the CMB power spectrum (and TE cross- 
correlation) . 

Each of our models is then characterised by the dark 
energy parameters Wde = {wo,Wm,at,A) and the cos- 
mological parameters Wc — {^DE,^bh^,h,ns,T,As), 
which are the dark energy density, the baryon density, 
the Hubble parameter, the scalar spectral index, the op- 
tical depth and the overall amplitude of the fiuctuations 
respectively. We therefore end up with ten parameters 
which can be varied independently. 

There is a remaining degeneracy in 715, r and flbh^, 
which allows the models to reach unphysically high values 
of the baryon density and the reionisation optical depth. 



Following the WMAP analysis and in order to remain 
consistent with |35j, we place a prior on the reionisation 
optical depth, r < 0.3. An alternative is to use a prior 
on flbh^, which may be physically better motivated and 
does a better job of breaking the degeneracy ~ we plan to 
use this for future work and discuss the difference later 
on. We also limit ourselves to models with w{z) > — 1, 
except where stated explicitly. 



III. EVALUATION OF THE LIKELIHOOD 

A grid-based method would necessarily lead to a very 
coarse sampling, which is why we opted for a Markov- 
Chain Monte Carlo (MCMC) method to sample the like- 
lihood surface. In addition this approach also allows easy 
marginalisation over parameters. We ran 16 to 32 inde- 
pendent chains on the UK national cosmology supercom- 
puter (COSMOS). This approach has both the advan- 
tage that there was no need to parallelise the Boltzmann 
solver, and that we are better able to assess the conver- 
gence and exploration by comparing the different chains. 

We take the total likelihood to be the product of the 
likelihoods of each data set (CMB, SN-Ia and LSS), or 
by defining xls = -21og£, 



Xtot 



XWMAP ■ 



■ XsNla(+X2dF)' 



(2) 



where the LSS contribution is added only in section lTV Dl 
We evaluate the WMAP likelihood using the code of the 
WMAP science team jSg , and we have checked that our 
SN-Ia likelihood results are consistent with ones of Riess 
et al |2J|. We treat the luminosity A^ as a nuisance 
parameter over which we marginalise analytically. This 
automatically also marginalises over the Hubble constant 
Hq-, so that the supernova data does not depend directly 
on it. For the 2dF results we use the formalism and data 
of Tegmark H^. 

The convergence and sampling behaviour of the chains 
in a high-dimensional space is far from trivial. Correla- 
tions in several dimensions are a particular issue as they 
lead to a high rejection rate. To improve the acceptance 
rate, we estimate in a first step the covariance matrix 
and then use rotated parameters which are linearly inde- 
pendent (i.e. lead to a diagonal covariance matrix with 
variances cr(pi)^ on the diagonal). We set the proposal 
width of parameter pi {i = 1, ..,N) to 2.5a{pi)/VN as 
advocated in |38ll39l|. which we found to work very well. 
We also evaluate the rejection rate every 100 steps and 
adjust the width dynamically, but normally this is not 
necessary once the covariance matrix is used. 

The result of a MCMC is a "cloud" of samples in each 
chain with a density proportional to the local value of 
the likelihood. In general we are not content with a N- 
dimensional description of the likelihood, but would pre- 
fer lower dimensional constraints. The generally accepted 
way to project the likelihood onto fewer dimensions is 
through marginalising. This means that we integrate the 



likelihood over the parameters which we want to hide, 
£{pi,...,Pi-i,p,+i,...,Pn) ex 
/ dpi£{pi, . . . ,pi^i,pi,pi+i, ...,p„) (3) 

This requires a measure fi — dp on the parameter 
space. In many cases, there is some physical motiva- 
tion for the choice of the measure. Alternatively, it may 
be that the for all reasonable choices, the measure does 
not vary strongly across the range of interest. As an ex- 
ample, we could use either fif, or flbh^ as a fundamental 
parameter. But since both Hq and flbh'^ are quite well 
constrained, the result will not change appreciably if we 
use the one or the other. In this case, it does not matter 
which one we use. If neither is true then the result can 
depend strongly on the choice of this measure. 

In the MCMC case, the choice of measure is implicit 
in the choice of parameters if one follows the usual rule 
that the marginalisation is done by summing up the sam- 
ples. To illustrate this, let us assume that the likelihood 
does not depend upon a given parameter 0.1 < p < 10. 
In this case, the marginalisation over this parameter is 
just the volume of the parameter space, since all values 
are equally likely. The resulting likelihood will be flat, 
independent of the choice of parametrisation. Using p 
as our fundamental parameter, we find that roughly half 
the points will be in p < 5 and half in p > 5. On the 
other hand, using logp we will find that half the points 
are in p < 1, and half in p > 1. The result depends in 
this case strongly on the choice of parametrisation. We 
would like to add here that a linear transformation of the 
parameters does not change the measure, as long as the 
boundaries of the integration are also adjusted. This is 
normally the case, since the integration volume is usually 
given by the region where £ 7^ 0. An example of such a 
transformation is the use of a covariance matrix to render 
the parameters linearly independent. 

A different way to think about this specific problem, is 
by looking at the choice of measure as a prior. Choosing 
for a parameter p the measure dp (i.e. sample it equally 
in p) is the same as imposing a flat prior, P{p) = 1. 
Instead, sampling the parameter in log(p), corresponding 
to a measure dp/p, corresponds to a flat prio r in log(p) 
or P{p) = 1/p when sampling evenly in p J40j . 

In order to test our dependence on the choice of 
parametrisation, we ran chains for different choices, 
namely (at. A), (log(at),A) and (log(at),log(A)). In the 
next section we will discuss which of the results depend 
on this choice. 

A second issue which is often neglected is that the like- 
lihood which we deal with here is by no means close to 
Gaussian in some of the variables. The structure of the 
minima can thus be arbitrarily complicated. It is there- 
fore important to specifically search for the global min- 
imum, which can of course be rather difficult in many 
dimensions. 

As a final remark before discussing the results, we 
would like to point out that caution is necessary when 



interpreting the results. One sigma limits are not suf- 
ficient to rule out anything - the true model is indeed 
expected to lie about one sigma away from the expecta- 
tion value! Even two sigma or 95% (the limits which we 
usually quote) are not sufficient. We present here limits 
and constraints on well over 20 different variables, and 
so we again expect at least one to lie in the excluded re- 
gion. Even worse, Gaussian statistics are well known to 
underestimate the tails of "real world" distributions, so 
that outliers are far more common than naively expected. 
In many fields, e.g. particle physics, a 5 sigma limit is 
being used to claim an actual detection. Although our 
data is not yet of sufficient quality to impose such strin- 
gent limits, we should bear this in mind. So we should 
as an example only consider regions ruled out in the 1- 
dimensional likelihood plots if the likelihood has (visu- 
ally) fallen to zero. 



IV. RESULTS 

In this section we discuss different aspects of our re- 
sults. We start by taking a look at the best-fit models and 
at the goodness-of-fit of both quintessence and ACDM 
models. Then we show that the introduction of a time- 
varying equation of state for the dark energy component 
does not significantly alter the constraints on the basic 
cosmological parameters. This allows us to discuss con- 
straints on the time evolution of the quintessence equa- 
tion of state. In section HVDI we use large-scale structure 
data instead of the supernova data to break the geomet- 
ric degeneracy. We also check if the combination of both 
data sets improves the constraints. Finally, we discuss 
limits on toy "phantom energy" models where we allow 
w < -I. 



A. The goodness of fit 

Our global best fit QCDM model is characterized by 
the following dark energy parameters: Wq — —1.00, 
Wm = —0.13, at = 0.48 and A — 0.06, which correspond 
to a fast transition in the equation of state at redshift of 
1. The total x^ oi this model is 1602.9, while the best- 
fit ACDM model has x^ = 1605.8. However, the total 
number of degrees of freedom is 1514, so that all our fits 
are rather bad. This is mainly due to the WMAP data 
(see the discussion in Q). In table U we report the corre- 
sponding x^ values for the CMB and SN data and best fit 
values of the standard cosmological parameters for these 
two models. Notice that the QCDM model provides the 
best fit to both the CMB and SN data. 

It is intriguing that such a model has a time evolving 
equation of state ^(z) similar to that reconstructed from 
the best fit to the SN data in [ij, IH E^ . In figure H 
we plot the temperature anisotropy power spectrum for 
these two models. It is remarkable how perfectly the 
two completely different models agree at intermediate 



Model XcMB XsN Xtot ™o ^de Hq Qth Us r 

ACDM 1428.7 177.1 1605.8 -1.0 0.69 69 0.023 0.97 0.11 
QCDM 1426.1 176.8 1602.9 -1.0 0.71 67 0.026 1.09 0.29 



TABLE L X Eind best fit values of the cosmological param- 
eters for ACDM and QCDM models. 



i, demonstrating the power of the WMAP CMB data. 
At low niultipoles the additional freedom of the QCDM 
models allows a slightly better fit. In fact due to a differ- 
ent distribution of the ISW, these QCDM best fit models 
have less power at low multipoles than the ACDM one. 
However we want to remark that this part of the CMB 
spectrum is most likely affected by galactic contamina- 
tion effects lU and without a more accurate investiga- 
tion we should not place too much emphasis on this sup- 
pression of power. It is worth noticing that the difference 
in the best-fit cosmological parameters between QCDM 
and ACDM will lead to different TT power spectra at 
higher multipoles, / > 700. This suggest that an accurate 
detection of the third peak may increase the statistical 
weight in favour or against the QCDM model. 
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FIG. 2: CMB power spectrum for the QCDM (red solid line) 
and ACDM (blue dashed line) best fit models. 



The fact that the x^ improves by 3 through the addi- 
tion of the three new dark energy parameters, is to be 
expected. Nevertheless we see that some quintessence 
models provide a better fit to the data than standard 
ACDM , as opposed to analyses which assume that w 
is constant. As we will discuss later, the limits on the 
time evolution of the dark energy equation of state are 
still compatible with a number of proposed scalar field 



scenarios. 



Studying the distribution of the x^ values in the 
MCMC chains for the ACDM models, we find that Ax^ = 
6.4 for the models at la (68.3% CL) and Ax^ = 11.8 at 
2cr (95.4% CL). In the Gaussian case, this corresponds 
to about 5.5 independent degrees of freedom, slightly less 
than the number of cosmological parameters used (6). 

For the quintessence models we find Ax^ = 9.9 for la 
and Ax^ = 15.5 for 2a. Assuming Gaussian errors, this 
would mean that we are dealing with about 8 indepen- 
dent degrees of freedom, many less than the 6 cosmolog- 
ical and 4 dark energy parameters used in the analysis. 
We conclude that we are unable to constrain all the ad- 
ditional parameters. 

Is the improved x^ of the best fit QCDM models a pos- 
itive evidence for the additional dark energy parameters? 
The answer to this question requires an estimation of the 
information criteria associated with the model |45l |. In 
general adding more parameters tends to improve the fit 
to data, however one should reward thoes models that 
can produce the same goodness-of-fit with fewer param- 
eters. This can be achieved through the Akaike informa- 
tion criterion (AIC l46l |) and Bayesian information crite- 
rion (BIG '47i|), respectively defined as 



and 



AIC 



BIC 



-2lnC + 2k, 



'2\nC + k\nN, 



(4) 



(5) 



where C is the maximum likelihood, k is the number of 
parameters of the model and N is the number of data- 
points. The AGDM models have an AIG of 1617.8, while 
the QGDM models of 1622.9. Based on this difference of 
5 we conclude that four parameter dark energy models 
are not necessarily favoured. The BIG disfavours the new 
parameters even more strongly. The information criteria 
therefore suggest that current data provide no positive 
evidence that the dark energy is anything more complex 
than a cosmological constant. This is not surprising, but 
should not be used as a reason to stop searching for ways 
of detecting evidence of evolution. The rewards from 
finding such evidence would be huge. 



B. Constraints on cosmological parameters 

The class of quintessence models could a priori contain 
new, severe degeneracies which change completely the 
preferred values of the cosmological parameters. If this 
were the case, then all the standard results of the WMAP 
analysis Ql would only be valid in the context of a AGDM 
model. We have found this not to be the case. The new 
dark energy parameters are degenerate amongst them- 
selves but do not introduce any new degeneracy with the 
other cosmological parameters. This can be seen in figure 
13 where quintessence model results are compared with 
those obtained for AGDM models only. Unless specifi- 
cally stated, the quoted results have been obtained using 



a linear parametrisation of at and a logarithmic one of 
A in the MGMG. This maximises the weight of models 
with a late rapid transition, other parametrisations lead 
to an even slightly better agreement between AGDM and 
QGDM. 
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FIG. 3: Marginalized likelihoods for the various cosmologi- 
cal parameters in the ACDM scenario (blue solid curve) and 
including the QCDM models (yellow shaded region). Also 
shown are QCDM models with a prior on the baryon energy 
density Qbh^ = 0.0216 ± 0.002 (red dashed line). The resuhs 
agree very well in all cases. 



We notice that all parameters, with the possible excep- 
tion of the reionisation optical depth r, are well deter- 
mined in both cases. Also, their values are very similar 
and cannot be distinguished even at one standard devia- 
tion. Table ITU lists the best fit values of the cosmological 
parameters for the AGDM and the dark energy models. 
The constraint on h is consistent with the HST measure- 
ment 30J , while the amount of clustering matter Q.^ is in 
agreement with large scale structure estimates [i^, Hy . 
The physical baryon density is consistent with the Big- 
Bang Nucleosynthesis expectations [3J|. This provides 
an important cross-check that all viable models have to 
pass. The background universe is therefore largely inde- 



Parameter 



ACDM 



QCDM 



n^ 0.30 ±0.03 0.29 ±0.04 

fifc/i^ 0.0237 ±0.0013 0.0240 ±0.0015 

Ho 70 ± 3 68 ± 3 

n^ 0.99 ±0.04 1.01 ±0.04 

T 0.15 ±0.06 0.19 ±0.07 



TABLE II: Mean and (formal) standard deviations for the 
cosmological parameters. The ACDM values agree with pub- 
lished analyses. The QCDM values are always consistent 
within one standard deviation, showing that quintessence 
does not significantly impact "precision cosmology" . 
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pendent of the details of the dark energy. 

We find nonetheless some differences, but they can be 
explained quite easily. Firstly, these are the marginalised 
likelihoods. The remaining degeneracy in fi™, Hq and wq 
is therefore translated into a slight shift to lower values in 
both rim and Hq. The degeneracy between the physical 
baryon density, flbh^-, and the scalar spectral index, ng, 
becomes slightly worse (see figure^, leading to some- 
what longer tails in their likelihood distributions. This is 
a consequence of the larger Integrated Sachs- Wolfe (ISW) 
effect produced in time dependent dark energy models 
with the respect to the ACDM case i28|. In fact the ISW 
boosts power on the large angular scales of the CMB, 
therefore less power of the primordial fluctuation power 
spectrum at small wavenumber k is required to match the 
data. Consequently slightly blue shifted spectral index 
values are preferred. The same eS'ect is also responsible 
for the spread of the likelihood distribution of the reioni- 
sation optical depth through the degeneracy between Ug 
and T. For larger values of Ug the excess of power at the 
low multipolcs of the TE spectrum requires larger values 
of the optical depth (see figure \^ . Since the reionisa- 
tion suppresses the contribution at high i with a factor 
of e"^"^ a better determination of the CMB peak struc- 
ture as well as better polarisation data will help to break 
this degeneracy. We hope that the latter will be available 
shortly, when WMAP releases the two-year data. 

Figure 121 shows also the effect of introducing a prior on 
flbh"^, namely Qth^ = 0.0216 ± 0.002. As we can see, the 
prior only affects the cosmological parameters r, ns and 
rifjft,^ and removes some of their high values. None of the 
dark energy parameters are affected (see figure EJ. We 
should note that all cases use a top-hat prior r < 0.3. 



C. Constraints on the dark energy parameters 

As we have seen in the previous section, the additional 
parameters which describe the dark energy do not in- 
troduce any new degeneracies with the standard cosmo- 
logical parameters. However we expect the dark energy 
parameter space to have an internal degeneracy. For in- 
stance Wo and Wm both act on physical observables as an 



FIG. 4: The degeneracy between the scalar spectral index, 
ns, and the physical baryon density Qbh^ ■ The filled contour 
are the 1 and 2 a limits of the quintessence models, while the 
black solid contours are those of the ACDM models. 




FIG. 5: The degeneracy between the scalar spectral index, 
ns, and the optical depth r. The filled contour are the 1 and 
2 a limits of the quintessence models. The black lines show 
the corresponding limits for the ACDM case. 



equation of state parameter. Therefore slowly varying 
models with Wm ~ wo niay look indistinguishable from 
models with the same value of wq and a rapid transition 
at very early time from whatever value of Wm- As antic- 
ipated in Section IIV Al the consequence of such degen- 
eracy is that we can strongly constrain only one of the 
dark energy parameters, which turns out to be wq. In 
particular, the one-dimensional marginalized likelihood 
gives Wo < —0.80 at 2a, which is consistent with the up- 
per limits quoted in other time dependent dark energy 
analyses [22, ISfl |51| ■ Of course what this really implies 
is that the equation of state at a redshift of order z=0.1 
is being constrained. We can not say anything about its 
true value today, but in what follows we use wq with this 



caveat in mind. The same applies when we take the Umit 
of z — > in the relevant figures. 

The other dark energy parameters are weakly con- 
strained. In particular, as expected from the arguments 
discussed in Section IIIII we find that the inferred con- 
fidence intervals may depend on the parametrisation of 
at and A used in the MCMC. For the standard case, we 
find Wm < —0.08 at 2a, on the other hand at and A re- 
main unconstrained. We refer to [23 for a discussion on 
breaking such a degeneracy with an estimate of the value 
of as from large scale structure. 

A parametrisation in log(at) would place more em- 
phasis on early transitions, so that the effective redshift 
where Wm is evaluated is moved to higher values. In this 
case the limits on Wm become even weaker, and we con- 
clude that Wm as a parameter is difficult to interpret. On 
the other hand, the limits on w{z) as a function of red- 
shift are less dependent of the parametrisation of at , as 
we derive them at all redshifts separately. We therefore 
advocate these limits, as plotted in figure as a better 
illustration of the constraints on dark energy models. 
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FIG. 6: Likelihoods for the dark energy parameters in the 
QCDM models without prior on Qbh^ (yellow shaded area) 
and with the prior (solid blue line). The prior does not affect 
the dark energy parameters significantly. We also show the 
relative goodness of fit of the best model in each bin (red 
dashed lines). 



These limits, derived with the Markov-Chain ap- 
proach, rest solely on the local density of the accepted 
models. If the assumption of Gaussian errors holds ap- 
proximately, then we can derive the same limits using the 
actual likelihood values instead. In this case, and consid- 
ering a single variable, models with Ax^ < 4 occur 95.4% 
of the time. We find that models with wq > —0.8 have 
Ax^ > 4, consistent with the limits from the Markov 



Chain. In order to put limits on the equation of state as 
a function of redshift, we proceed as follows: we compute 
w{z) for each model and then take either the 95% con- 
fidence region over all models in the chain (shaded area 
in figured) or compute the highest w{z) for models with 
X^ < Xmin + 4 (dashed line). 




FIG. 7: 2a confidence region on w{z) derived by taking the 
95% models with lowest w{z) from our main chain (yellow 
shaded area) and the "exclusion zone" where all models have 
Ax^ > 4 from the best-fit model in our main chain (blue 
hatched area). We also show the 'w{z) of the best-fit model 
(black solid line). ACDM is acceptable at 2a. 



The shaded area is a priori the proper marginalised 
result, representing the 95% limits on wiz). However it 
is worth remarking that this constraint may suffer from a 
potential problem with the choice of the measure that we 
introduce to integrate the marginalised parameters over 
(see discussion in section llTT|) . On the other hand the 
dashed line relies only on the goodness of fit, and can be 
interpreted as an exclusion limit. In other words mod- 
els with a w(z) that enter the hatched area are a "bad" 
fit to the data. This does not include any information 
about how likely these models are, given the variation 
in the other parameters. Our interpretation of this is as 
follows: if we want to judge if a single, specific model 
is ruled out or not, then we should be using the dashed 
line as an upper limit for w{z). If we are more interested 
in what we would expect as the value of w{z), given our 
additional knowledge about the other variables, then we 
should look at the shaded region as providing a limit on 
the equation of state parameter. As we can see in figure 
13 models with Wm > at z > 1 and with a fast transition 
occurring at z < 1 are a bad fit to the data, being beyond 
the dashed blue line. Physically this is because models 
where the transition from Wm to wq occurs at redshifts 



z < 10 with Wm > —0.1 give rise to a non-negligible 
dark energy energy contribution at decoupling, which is 
strongly constrained by CMB data. However, models 
with wq < —0.8 and Wm > —0.1 for which the transi- 
tion occurs at redshift z > 10 are consistent with the 
data as their early energy contribution is negligible. As 
mentioned earlier, the apparent exclusion of these models 
based on the likelihood for Wm in figure |3 is an artifact 
of our parametrisation of at- This can also be seen by 
noticing that the maximised results (dashed curve) does 
not fall to zero for w^ > 0, indicating that there are 
acceptable models in this part of the parameter space. 
However these models will be indistinguishable from a 
pure ACDM scenario and are thus not very interesting 
when we try to rule out one from the other. 

The limits on w(z) can be reinterpreted as constraints 
on the evolution of the dark energy density (see figure|Hl • 
Models with flnEiz) above the dashed line are ruled out, 
thus limiting the amount of dark energy available during 
matter domination to be ^Dsiz) < 0.1 (see also [SOllS^ '). 

These results undoubtedly have an effect for 
quintessence model building, and we will be investigating 
this aspect in a future paper. 
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FIG. 8: Limits on Q,de[z) corresponding to those on w{z) in 
figure 



For now though we restrict ourselves to a few gen- 
eral remarks. The limits on w{z) previously obtained 
allow us to constrain a large class of quintessence mod- 
els. For instance the exclusion plot in figure suggests 
that models with a perfect tracking behavior, for which 
w = during the matter era up to z < 4, and with a late 
time fast transition, are disfavoured by the data. As we 
have seen before, this is because fig must be negligible 
at early times. This class of models, for some partic- 
ular choice of the parameters of the scalar field poten- 



tial, include the two exponential potential |53| and the 
Albrecht-Skordis model |5J| . Of course, if they leave their 
tracking behaviour before then, the constraint is weak- 
ened. These models satisfy the constraints only if the 
slope of the scalar field potential, where the quintessence 
is initially rolling down, is very steep and then followed 
by a nearly flat region such that the equation of state 
sa — 1 at the present time. On the other hand models 
with a non-perfect tracking behavior and a slowly vary- 
ing equation of state with wq < —0.8 are consistent with 
the data. This is the case of quintessence models with an 
inverse power law potential [3] , super gray ity inspired po- 
tentials such as the one proposed in [53 or off-tracking 
quintessence models, such as those studied in |^ l57l |. 
Models of late time transition ^«| which show features 
of our best fit model can also be consistent with the data. 



D. Adding large scale structure 



One important question is whether the supernova data 
contains any severe systematic effects that may strongly 
bias our results. To test for this possibility, we replace 
the supernova data by the 2dF galaxy redshift survey 
power spectrum from ref. [37| . The bias is added as a 
new, free parameter. This means that only the form of 
P{k) is constrained, not its amplitude. The solid lines in 
figuresl^andllOlshow the resulting likelihoods of the cos- 
mological and the dark energy parameters respectively. 
Comparing them to the CMB-|-SN-Ia results (shaded yel- 
low regions), we find that the results are consistent, but 
that the supernovae allow us to place stronger constraints 
on r^m, Hq and wq when combined with the CMB. There 
is no evidence for any systematic problems in the SN-Ia 
data set. 

Additionally, we can use the combination of all three 
data sets. In |23| we found that the large scale structure 
(LSS) does not add any strong constraints, beyond those 
found with CMB-|-SN-Ia data, as long as no constraints 
on the bias parameter are imposed. As the dashed curve 
in figures and ^] shows, this is still the case, and the 
bias parameter is strongly correlated with the cluste ring 
strength, a^. The constraints on the bias found in [53 
do not apply to our analysis, since they were obtained by 
combining the SDSS 3D matter power spectrum with the 
WMAP results on erg which are correct only for ACDM 
cosmologies. 

For standard quintessence models one can thus either 
use supernovae or large scale structure data, with the 
supernova data giving the stronger constraints. The sit- 
uation changes if additional parameters need to be con- 
strained, e.g. non-zero neutrino masses. In this case it is 
crucial to have very good clustering data on small scales 
where the neutrinos impose a distinct signal on P{k). 
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FIG. 9: Constraints on the cosmological parameters when 
adding the 2dFGRS data to CMB and SN-Ia data (blue 
dashed line) and when only using CMB+2dFGRS (black solid 
line), compared to the CMB+SN-Ia case (yellow shaded re- 
gion). 



E. Phantom energy models 

Several authors have suggested that dark energy mod- 
els with a super-negative equation of state, for which 
w < —1, can provide a better fit to the CMB and the SN- 
Ia data jMimileliJ. On the other hand all these 
analyses are biased in favour of phantom dark energy 
models since they use a constant equation of state pa- 
rameter 1(34, "65]. In this scenario the dark energy violates 
the weak energy condition (WEC) which leads to a num- 
ber of problems JBg . For this reason we feel that these 
models are disfavoured on theoretical grounds. Neverthe- 
less it remains an interesting question whether they are 
compatible with the current cosmological observations? 
Constraining time dependent phantom energy models al- 
lows us the opportunity to test these models against the 
observational data without the bias induced by assum- 
ing a constant equation of state parameter. However the 
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FIG. 10: Constraints on the dark energy parameters when 
adding the 2dFGRS data to CMB and SN-Ia data (blue 
dashed line) and when only using CMB-f 2dFGRS (black solid 
line), compared to the CMB-|-SN-Ia case (yellow shaded re- 
gion). 



main problem for this class of models is the existence of 
tachyon instabilities which lead to an exponential growth 
of phantom energy perturbations on small scales. In 
addition our formalism does not allow us to follow the 
evolution of the dark energy perturbations for phantom 
models which cross the w = — 1 value (see discussion in 
Appendix IB|| . On the other hand we can account for the 
perturbations in models which violate the weak energy 
condition at all times. 

In order to be self-consistent we therefore extend our 
analysis to two different classes of dark energy models, 
those which always satisfy/violate the WEC and those 
which cross the WEC boundary value w = —1. The lat- 
ter are assumed to be homogeneous and consequently our 
analysis for these class of models accounts only for the 
effects they produce on the background expansion. Al- 
though this is not physical, we are unaware of a unique 
prescription for handling these cases. We therefore sug- 
gest that the reader sees this section more as a specu- 
lative treatment. It is interesting to note that in these 
"toy" phantom models, the allowed values of the cosmo- 
logical parameters do not change very much, in that they 
lie somewhere between the ACDM and the QCDM case 
(seeFig.inj- 

The models which cross w — —1 provide a slightly im- 
proved best-fit. In particular we find a best-fit model 
with x^ — 1601. It has wq = —2.0 and Wm = —0.7, 
i.e. the equation of state crosses over that of the cosmo- 
logical constant, p — —p. This behaviour is most likely 
driven by the supernova data, as we find a similar result 
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FIG. 12: Constraints on the dark energy parameters for all 
phantom models (red dashed lines) and those that respect 
the weak energy condition (w > — 1) (solid black line). The 
dotted line is from models that do not cross w = —1. 



regions correspond to all models which do not cross w — 
— 1 and for which perturbations are included, while the 
hatched area corresponds to all models which cross w = 
-1. 



FIG. 11: Constraints on the cosmological parameters for all 
phantom models (solid lines) and those that respect the weak 
energy condition (w > 1) (shaded), compared to the ACDM 
models (dashed). 



if we only use this data set |27| . 

For the standard parametrisation of at and A, the me- 
dian for wq is —1.1 and for Wm is —0.8. The 95% con- 
fidence intervals are —2.7 < wq < —0.77 and —8.3 < 
Wm < —0.19. For the models that do not cross w = —1, 
but always remain either with «;<- lorw;>— 1, the re- 
sults give the median for wo to be —1.0 and for Wm to be 
—0.97. The overall best-fit model is the same as for the 
standard QCDM models, but as the above median val- 
ues show there are about the same number of accepted 
models on both sides of the divide. 

Thus there is no evidence for any deviation from 
ACDM in this extended framework. 

The parameter which does show some change is the 
clustering amplitude, as, which lies now in the 95% con- 
fidence interval 0.66 < era < 1-24 as opposed to the 
dark energy models satisfying the weak energy condition 
w> -1 for which 0.53 < cts < 1-07. 

The limits on w{z) are shown in figure [T^ The shaded 



V. CONCLUSIONS AND OUTLOOK 

In this paper we have made the first attempt to con- 
strain dynamical models of dark energy by combining 
CMB, 2dFGRS and type-la supernova data. On the pos- 
itive side, we find that by allowing the dark energy equa- 
tion of state to vary as a function of time, we do not intro- 
duce any strong degeneracies that would adversely affect 
the standard cosmological parameter estimation, apart 
from perhaps a mild degeneracy between the reionisation 
optical depth, r and the redshift of the commencement of 
acceleration. Its effect is to alter the ISW contribution to 
large angles and hence, after COBE normalisation uni- 
formly alters the heights of the peaks Q2^ . Breaking this 
degeneracy by some means will significantly enhance our 
ability to constrain dark energy dynamics. This might 
be done through a number of routes: other astrophys- 
ical constraints on r; by probing the redshift of accel- 
eration using ISW-LSS correlations; by measuring the 
non-Gaussianity of CMB from weak lensing induced by 
structure formation [63| or by beating cosmic variance 
using cluster polarisation at high redshift |63| ■ 

The remaining cosmic parameters are only mildly af- 
fected by the new freedom given to dark energy and are 
similar to their ACDM counterparts. On the down side. 
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unstable, and so we turned them off for any model which 
enters w < —1 at any point. In this case models that have 
w slightly larger than —1 at early times and show a rapid, 

late transition to a super-negative equation of state, w '^ 
—2, are slightly, but not significantly, preferred. In a 
slightly more physical model where w = —1 cannot be 
crossed, allowing us to compute the perturbations, we 
find no preference for these "phantom" models. 

Overall, we conclude that dynamical models of dark 
energy are by no means ruled out and provide a slightly 
better fit to current data than ACDM . The latter is 
still perfectly acceptable and, given its simplicity, seems 
in many ways the preferred model at the current time. 
Whether observations and theoretical prejudice maintain 
this conclusio over the coming years will have a profound 
impact on our understanding of fundamental physics. 
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FIG. 13: Constraints on the dark energy equation of state pa- 
rameter w{z). The shaded area corresponds to models which 
do not cross w = — 1 and for which perturbations are taken 
into account, while the hatched area corresponds to all mod- 
els, including those crossing w = —1, but without perturba- 
tions. 



with regard to the dark energy parameters, the currently 
available data set means that only the present value of 
the equation of state wq is well constrained. We find 
Wo < —0.8 at 95% confidence level. By studying the be- 
haviour oiw{z), we conclude that the constraints become 
weak for any redshifts larger than about unity. 

However, there are some clear results that emerge. For 
example in the large class of Quintessence models that 
have periods of perfect tracking behaviour, i.e. w = 
during the matter dominated era, only those which track 
for during the period z ~ 5 are acceptable. Quintessence 
models with a very late departure from tracking, or any 
dark energy models with a late transition from a high 
value of w, are disfavoured as a consequence of the fact 
that early contributions of dark energy are constrained 
to be neghgible. 

We have also found a point of practical importance in 
the hunt for dark energy dynamics. Since only rapidly 
varying models at low redshift have a distinct signature 
as opposed to ACDM [23 including such models in ones 
likelihood analysis is very important, otherwise the re- 
sults will be biased towards no detection of dynamics. 
However these rapidly varying models are the most sus- 
ceptible to numerical errors, which we have found can 
be as large as 5%. We therefore have used a code spe- 
cially adapted to handle these "kink" cases, and this is 
described in Appendix [HI 

We have also studied toy models for which w < — 1 is 
possible. In this case the perturbations are generically 
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APPENDIX A: PARAMETRISATION OF THE 
MODELS 

A simple way of describing a dark energy component 
is to consider a constant equation of state w. However 
such an approach suffers from a number of drawbacks. In 
fact assuming w to be constant introduces a bias in the 
analysis of cosmological distance measurements with the 
consequence that large negative values of w are favoured 
if the dark energy is time dependent [63, [63 • The effect 
of this bias has to be carefully taken into account partic- 
ularly when models of phantom energy |6(ll66j|. for which 
w < ~1, are constrained. A possible way to address this 
issue is to constrain a time dependent parametrisation 
of w{z). In fact the time dependence of the dark en- 
ergy equation of state completely specifics the evolution 
of the dark energy density through the continuity equa- 
tion. Several formula have been proposed in the litera- 
ture [S IZS El El all with limited applicability. They 
are typically based on Taylor expansions in some variable 
(e.g. z, log(l -I- z) or 1 — a). 
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However, tracking quintessence models exhibit the im- 
portant property that before the universe begins to ac- 
celerate the dark energy has an equation of state which 
mimics that of the dominant energy component. Today 
of course w < —1/3, so there must be a transition be- 
tween the value at high redshifts and that at the present 
time at some critical redshift, zt, or equivalently scale fac- 
tor, at = 1/{1 + Zt), with a thickness determined by a pa- 
rameter A. In fact, this physically motivated parametri- 
sation bears a striking resemblance to the tanh(z) kink 
soliton solution in spirit. 

To compare various proposals we write it as: 



discussion). For instance, models characterized by a 
slowly varying equation of state, such as supergravity 
inspired models [53, correspond to a region of our pa- 
rameter space for which < at/A < 1, while models 
with a rapid variation of w(a), such as the two expo- 
nential potential [53 or the Albrecht-Skordis model |54| . 
correspond to at/ A > 1. Models with a simple constant 
equation of state are given by wq = Wm- We can ac- 
count also for the so called phantom energy models for 
which WQ,Wm < —1. The cosmological constant model 
corresponds to the following cases: wq = Wm = — 1 or 
Wq — —1 and at ^ 0.1 with at/A > 1. 



w{a) = Wo + (wm - wa)V{a, at, A); 



(Al) 



where F determines precisely how the transition from 
w = Wm to w = Wq occurs. Two proposals for F have 
been made in the literature. The original one |17ll2^ uses 
z instead of a: 



T{z,zt,A) 



-1 



1 + e- 



(A2) 



while a different proposal (and the one that we use here) 
was made in |26J, namely: 



F(a,at, A) 



1 + e* 
1 + e-^ 



1 



1 



eA 



(A3) 



The latter form has the advantage of being more stable 
numerically (since a is bounded in the interval [0, 1] while 
redshifts up to z = 10^ must be considered in the first 
form) and satisfies w{a = I) = wq whereas the same is 
only true for sufficiently rapid transitions in eq. I|A2|I . 
Because of this we use the latter form (|A3|) in this paper. 

Importantly, this form can also be extended 26] to al- 
low for a different value of w during radiation domination 
so that it can represent tracking models faithfully at all 
redshifts in terms of five physical parameters: the value 
of the equation of state today wq, during the matter era 
Wm and during the radiation era w,-, while the time de- 
pendence is specified by the value of the scale factor at 
where the equation of state changes from Wm to wq and 
the width of the transition A. Since Big-Bang Nucle- 
osynthesis bounds limit the amount of dark e nerg y to be 
negligible during the radiation dominated era [73 the ex- 
tra freedom is not particularly important and we further 
reduce our parameter space by setting Wr = Wm yielding 
the transition function HA3|I . 

Written out fully then, the equation of state we use in 
this paper is Eq. (4) of ref . _26J : 



w{a) = wo + iwm~wo)x 



1 



1 



1 + e- 



1 



(A4) 



The dark energy parameters specified by the vector 
W DE = {wQ,Wm,a't,A) with the parametrisation given 
by Eq. (|A4|) can account for most of the dark energy 
models proposed in the literature (see [23 for a detailed 



APPENDIX B: COSMOLOGICAL EVOLUTION 
OF SCALAR FIELDS 

The cosmological evolution of minimally coupled 
quintessence field Q is described by the Klein-Gordon 
equation. 



+ 2-Q' 






0, 



(Bl) 



the prime denotes derivatives with the respect to confor- 
mal time, V{Q) being the scalar field potential and 



SttG 



2a2 



V{Q) 



(B2) 



where pm and pr are the matter and radiation en- 
ergy density respectively. The equation of motion for 
the quintessence fluctuations at the scale k in the syn- 
chronous gauge is given by 



SQ" + 2-5Q' + {a^^ 
a dQ'^ 



1 



k^)5Q = --Q'h', (B3) 



where h is the metric perturbation. Instead of specifying 
the scalar field potential V{Q) the conservation of the 
energy momentum tensor allows us to describe a scalar 
field as a perfect fluid with a time dependent equation 
of state w{a). In such a case the dark energy density 
evolves according to 



PDE{a) 



3Hl 

SvrG 



Ude exp 



1 -f- w{s) 



ds 



(B4) 



where Hq is the value of the Hubble parameter and ^de 
is the dark energy density. We can describe the pertur- 
bations in a dark energy fluid specified by w{a) using 
Eq. IJB3|) with the second derivative of the scalar field 
potential is written in terms of the time derivatives of 
the equation of state w{a) |73| . 



,fv_ 



-|<^ 



a 
a 



7 3 

2 + 2^ 



w 



w ,a' 

Ail + w) 2 a^ ' 



(B5) 
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FIG. 14: Relative difFerence of CMB spectrum for a 
rapidly varying dark energy model computed with CM- 
BEASY, KINKFAST 1.0.0 and CMBFAST 4.5.1. 



For models with an equation of state rapidly evolving to- 
wards w — —1 the Eq. (|iJ5|l is undefined since w = —1 
and w'jw" = after the transition to —1. This corre- 
sponds to the fact that at the classical level the vacuum 
state has no perturbations and quantities such as the 
sound speed are not defined anymore. When such con- 
ditions are realized by one of our models in the MCMC 
chain, we set the dark energy density perturbations to 
zero by the time at which w — —I. On the other hand 



Eq. (|B5|I becomes singular at w = —1 if w' , w" ^ 0, which 
could be the case for phantom dark energy models with 
w{a) crossing the value —1. On the contrary perturba- 
tion in a time dependent phantom fiuid with w < — 1 at 
all the time can be described in our framework by switch- 
ing the negative sign in front of the second derivative of 
the potential in Eq. ljB3p . 



APPENDIX C: THE BOLTZMANN CODE - 
KINKFAST 1.0.0 

We have modified the CMBFAST 4.5.1 to implement 
our parameterisation of dark energy (KINKFAST 1.0.0). 
We have tested the numerical accuracy of our code by 
comparing the power spectra of different dark energy 
models with those computed by using CMBEASY |75l| . 
We find a perfect agreement within 1%. On the other 
hand we notice that the current version of CMBFAST 
4.5.1, which implements the dark energy by reading a 
redshift sampled equation of state, leads to wrong spec- 
tra in the case of dark energy models with a rapid evo- 
lution in the equation of state (see figure EJ. The cause 
of such a discrepancy is that the dark energy density is 
obtained through a splint integration procedure of the 
sampled equation of state. We find this method to be 
of poor accuracy. Without an analytic formula for w{z), 
it is more reliable to derive a polynomial fitting func- 
tion for the sampled equation of state and integrate the 
corresponding polynomial form with standard Numerical 
Recipes |76| integration subroutines such as rombint 
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